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The nature of the dark energy is still a mystery and several models have been proposed to explain 
it. Here we consider a phenomenological model for dark energy decay into photons and particles as 
proposed by Lima [1]. He studied the thermodynamic aspects of decaying dark energy models in 
particular in the case of a continuous photon creation and/or disruption. Following his approach, 
we derive a temperature redshift relation for the CMB which depends on the effective equation of 
state We// and on the "adiabatic index" 7. Comparing our relation with the data on the CMB 
temperature as a function of the redshift obtained from Sunyaev-Zel'dovich observations and at 
higher redshift from quasar absorption line spectra, we find w^ff = —0.97 ± 0.034, adopting for 
the adiabatic index 7 = 4/3, in good agreement with current estimates and still compatible with 
We// = — 1, implying that the dark energy content being constant in time. 
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■j-j , After the discovery that the cosmic expansion is accelerating [2, 3] and the first cosmic microwave background 

^ ' radiation (CMB) observations of a flat universe [4], the current standard model of cosmology implies the existence 
of dark energy which accounts for about 70% of the total energetic content of the universe, which according to the 
observations is spatially flat [5]. The nature of the dark energy is still a mystery (see for instance the review of [6]). 
Several models have been proposed to explain dark energy [7-15]. An alternative consists to consider a phenomeno- 
logical decaying dark energy density with continuous creation of photons [16, 18] or matter [15]. The dark energy 
^1^1 might decay slowly in the course of the cosmic evolution and thus provide the source term for matter and radiation. 

fT^ , Different such models have been discussed and strong constraints come from very accurate measurements of the CMB. 
P\| ' CMB is the best evidence for an expanding Universe starting from an initial high density state. Within the 

r [ Friedmann-Robertson- Walker (FRW) models of the Universe the radiation after decoupling expands adiabatically 

f~^ and scales as (1 + z), z being the redshift. Depending on the decay mechanism of the dark energy the created photons 

f^ ' could lead to distortions in the Planck spectrum of the CMB. If the dark energy is considered as a second fluid 

^-H . component transferring energy continuously to the material component, the second law of thermodynamics constrains 

►^ ' the whole process in such a way that the temperature law, that is how the temperature of the CMB changes with 

time, can be determined. Such a model has been discussed in detail by Lima [1], who also established under which 

conditions the equilibrium relations are preserved. He also introduced the concept of an "adiabatic" vacuum decay. 

'^ In a further paper Lima et al. [16] following the approach outlined in [1] derived a temperature redshift relation law 

■ - T ~ {1 + zY~^ assuming a nonvanishing source term in the balance equation for particle number, which depends on 

a function j3 assumed to be constant. This parameter can then be determined by fitting the temperature redshift 

relation to the data. However, being a purely phenomenological parameter it cannot be related to more fundamental 

quantities such as for instance the effective equation of state of the dark fluid Weff ■ 

In this paper we follow the same approach by Lima [1] , we do not restrict to the case of a universe dominated by 
a radiation fluid as in [16], since we leave the adiabatic index 7 (which appears in the equation of state of the fluid 
including both matter and radiation) free to vary and in particular we do not adopt their simplifying assumption 
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of a constant /3 and thus we derive a temperature redshift relation for CMB which depends on more fundamental 
quantities such as Weff and 7. We compare our temperature redshift relation with the available data on the CMB 
temperature as a function of the redshift obtained from Sunyaev-Zel'dovich observations and at higher redshift from 
quasar absorption line spectra. 

The paper is organized as follows. In § II wc will introduce the theoretical background for our model, deriving the 
main quantities of interest. In § III the multi-rcdshift measurements of the CMB temperature are presented, while 
§IV and §V are devoted to the discussion of the results, conclusions and future prospects. 

II. TEMPERATURE-REDSHIFT RELATION 

We assume the usual Friedmann- Robertson- Walker (FRW) cosmology and that the material medium of the Universe 
contains three different components: a fluid, which includes radiation and matter (both baryonic and dark matter) 
as particular cases, for which we have the law 

p={j-l)p (1) 

with 7 in the interval [0 - 2]. Further we assume a dark energy, quintessence-like x component, with pressure p^ and 
density p^ and a 'bare' cosmological constant Aq. With these components we get for the Einstein field equations 

87TGip + p,)+Ao = 3^+3^ , (2) 

87rG(p + p.)-Ao = -2|-|J-A, (3) 

where a dot means time derivative. For the dark energy we assume the relation p^ = —p^ (thus setting Wx = —1) to 
hold. This is for instance the case for quintessence models in the limit where the scalar field does not depend on time 
and thus its time derivative vanishes. In the later stages of the Universe the time dependence is possibly very weak so 
that Wx = —1 holds up to small corrections, which we will neglect in the following, thus simplifying our calculations. 
This can be viewed as a first order approximation. 

Furthermore, we assume that there is no curvature, thus k — Q, and take the sum of eqs. (2) and (3). This way we 
get 

8^G(p + p) = 2|^-2| = -2iJ. (4) 

Following the paper by Lima et al. [16] the energy conservation equation can be written as 

p + ?,{p + p)H = Cx, (5) 



where H — R/R is the Hubble parameter and 



Cx = -px -HPx+Px)H , (6) 

is a term which depends on the dark fluid and acts as a source term for the 7-fiuid energy. Evidently, if no interaction 
between the different fiuids exists, then Cx is null and the standard picture is recovered. Cx can describe different 
physical situations such as for instance a thermogravitational quantum creation theory [17] or a quintessence scalar 
field cosmology [13]. 

Assuming as mentioned above the relation px = —px and writing px ~ A(t) / (SnG) , we have 

On the other hand the equation for the particle number density is given by 



h + 3n_ff = i/i 



(8) 



where n is the particle number density and tp is the particle source term. Using Gibbs law and well-known thermo- 
dynamic identities, following the derivation given in the paper by Lima et al. [16], one gets 
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To get a black-body spectrum the second term in brackets in eq.(9) has to vanish, thus 



n 



(10) 



Note that Cx is null if p = — p, (this is the trivial case when only a dark fluid exists) or if i/i = 0, which means that 
both the dark and normal fluids are separately conserved. This way we get with eq.(l) 



and eq.(9) becomes 
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With (^\ = (7 - 1) one obtains 
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Using the equation for the particle number conservation eq.(8) into eq.(13) leads to 
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With eqs.(lO) and (7) we get 
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We insert (4) in eq. (15) and obtain 
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As next we integrate cq.(16) 
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where ig denotes the present time and ti some far instant in the past. With H = R/R the second term in the bracket 
on the right and the term on left hand side can be immediately integrated. Indeed, if A vanishes and 7 = 4/3 one 



gets the usual dependence T{t) — ^ ^j^/^!^ ''' for a radiation fluid. Whereas to carry out the integration of the first 
term on the right hand side it is useful to perform a change of variable from t to z and accordingly ^ = h(i+ ) ■ This 
way we get (with zi corresponding to the time ti and zq = corresponding to to present time) 
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where ' denotes derivative with respect to z. We now assume a power law model for the A term, thus A ~ B(R/Ro)^™ 
or A = B{1 + z)™, where i? is a constant, which is i? = 3i?o(l ~ ^mo) [15], if the "bare" cosmological constant Aq 
vanishes. 

For such a model the Hubble parameter as a function of m can be computed [15] and leads to 



H{z) = Ho 



(3r^mO - m) 
3 — m 
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where flmo is related to the present value of the total matter density: pmo ~ ^H^Mpiflmo, with Mpi — [SttGY 
Due to the negligible energy density accounted for by the radiation, this term has been neglected here. As next we 
insert H{z) and its derivative as taken from eq.(19) into eq.(18) and integrate it, to get (setting zi ~ z) 
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where 



A = [ln{{m - 3a„o) + ™(1 + z)"-^(f^™o - 1)) - ln{{m - 3)n„o)] 
We can also write eq.(20) as 
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We inserted in the exponent of eq.(20) the explicit form of B, thus getting as exponent in the above eq. (7 — 1). 
Notice that for z = 0: T(0) = To, whereas for to = the expression in the parenthesis is equal to 1 and thus 
r(z) = To(l + z)^^'*~^\ which for the canonical value of 7 = 4/3 reduces to the standard expression. Also for to, = 3, 

although not a realistic value, the expression in the parenthesis is non-singular and equal to: 1 H — — ^^ "*" — -. 

Since from the analysis by Ma [15] we expect m to be small we give here an expression for T(z), where we keep 
only the leading linear terms in the expansion in to and neglect all terms of higher order. This way we get from eq. 
(22) 



r(z) = To(i + z)^(^-i) (1 + y - ^(1 + (1 + ^r'i^mo 1))) 



(7-1) 
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If m is positive, then the dark energy slowly decreases as a function of the cosmic time, whereas if m is negative 
the inverse process happens. 

It is easy to calculate the equivalent effective dark energy equation of state p = w^ffP with w^ff = y — 1. If7n>0 
then we have Weff > —1, i.e. our model is quintessence-like [9, 12, 13], while we have a phantom-like [11] model when 
TO, is negative and iCe// < — 1- Another interesting quantity is the deceleration parameter, which can be written as 
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(24) 



Imposing that q{z) = 0, we can determine the transition redshift, i.e. the redshift when the Universe was changing 
from a deceleration to an acceleration phase, which is given by 
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- 1 



3flmO - m 
From this result we have that the larger m is, the earlier the Universe changes from deceleration to acceleration. 



(25) 



III. MULTI-REDSHIFT MEASUREMENTS OF Tcmb AND FITTING PROCEDURE 

To test the decaying A model we will rely on the CMB temperatures derived from the absorption lines of high 
redshift systems and the ones from Sunyaev-Zel'dovich effect in clusters of galaxies (we will collectively quote as 
Tcmb, hereafter). 

We can infer the CMB temperature at high redshift from the analysis of quasar absorption line spectra which give 
atomic or ionic fine structure levels excited by the photo-absorption of the cosmic microwave background radiation. 
Detection of absorption from the ground states and excited of Ci in damped Lya system towards quasars permits 
to estimate the population ratio of the excited fine-structure levels, and thus to derive the CMB temperature at the 
redshift of the damped Lya system. Ge et al. [19], in the damped Lya system of the QSO 0013-004, found the CMB 
temperature 

Tcmb = 7.9 ± 1.0 K at z = 1.9731. (26) 

Srianand et al. [20] detected absorption lines in an isolated gas cloud toward PKS1232-)-0815, with which they 
could get the following limits on the CMB temperature 

6.0 < Tcmb < 14 K at z = 2.33771, (27) 

and Molaro et al. [21] in the damped Lya system toward QSO 0347-3819 gave 

Tcmb = 12.113:2 K at z = 3.025. (28) 

The cosmic microwave background can also excite levels of molecular species, when the energy separation involved 
corresponds to the CMB peak frequency. We know that molecular gas is an important ingredient of star formation. 
Damped Lya absorbers (i.e. DLAs) are generally taken as seeds of present-day galaxies, and for this reason DLAs 
can be considered as an important gas reservoir for star formation. H2, the most abundant molecules in the Universe 
is a good candidate (see Puy et al. [22], Galli & Palla [23] and Standi et al. [24]). The search of H2 in DLAs can 
be carried out by observing the H2 absorption lines in the lyman X^S+ 1 — > B^ S+. Cui et al. [25] obtained the 
spectrum of molecular hydrogen associated with the damped Lya system at z = 1.7765 toward the quasar Q1331-I-170 
and constructed a model to describe the structure of H2 absorber. Applying the inferred conditions to the Ci fine 
structure excitation, they found the CMB temperature to be 

Tcmb = 7.2 ± 0.8 K at z = 1.7765. (29) 

Srianand et al. [26], from the CO rotational excitation temperatures, derive in a damped Lya system towards SDSS 
J143912. 04+111740.5 

Tcmb = 9.15 ±0.72 K at z== 2.41837. (30) 

During passage through a cluster of galaxies some of the photons of the cosmic microwave background radiation are 
scattered by electrons in the hot intracluster medium. This imprint was first described by Sunyaev-Zel'dovich [27]. 
Thus, spectral measurements of galaxy clusters at different frequency bands yield independent intensity ratios for each 
cluster. The combinations of these measured ratios permit to extract the cosmic microwave background radiation 
(see Fabbri et al. [28]). Recently, Luzzi et al. [29] have analyzed the results of multifrequency Sunyaev-Zel'dovich 
measurements toward several clusters from 5 telescopes (BIMA, OVRO, SUZI II, SCUBA and MITO), see Table 1. 

We will match the observed Tcmb with the theoretical expression Tj/j, which we have derived in eq. (22), by 
minimizing the following merit function 

N tCM B / rpi _ rpi \ 2 

XICMB= E *" ^^^ ) ' (31) 

■^ \ CFCMB.i / 

where (JcMB,i is the error on the temperature estimates and Ntcmb ~ 18 is the number of available observational 
data. We will find the best fitted parameters which correspond to the minimum of x^, i-e. Xmin ^^^ '^^ will determine 
the 68% uncertainties by imposing Ax^ — X^ — Xmin = 1- 



TABLE I: CMB Tcmb values for several clusters from Luzzi et al. [29]. 



Cluster 


Tcmb 
(K) 


z 


A1656 


2.72 ± 0.10 


0.023 


A2204 


2.90 ± 0.17 


0.152 


A 1689 


2.95 ± 0.27 


0.183 


A520 


2.74 ± 0.28 


0.200 


A2163 


3.36 ± 0.20 


0.202 


A773 


3.85 ± 0.64 


0.216 


A2390 


3.51 ± 0.25 


0.232 


A1835 


3.39 ± 0.26 


0.252 


A697 


3.22 ± 0.26 


0.282 


ZW3146 


4.05 ± 0.66 


0.291 


RXJ1347 


3.97 ± 0.19 


0.451 


CL0016+16 3.69 ± 0.37 


0.546 


MS0451 


4.59 ± 0.36 


0.550 



IV. RESULTS 

We have tested our model by comparing the CMB temperature predicted (see eq. 22), vi^ith the collection of multi- 
redshift measurements of Tcmb we have discussed in the previous section. We set Tq — 2.725 K, which is quite well 
determined in the literature [30], and the matter density fimo = 0.273 to the value inferred in Komatsu et al. [31]. 

If we take 7 = 4/3, then we find m — 0.09 ± 0.10 and zt = 0.82 ± 0.10, while, if we leave both 7 and m free to 
vary m = 0.201q23 and 7 = 1.35^0 05 and zt ~ 1-02 ± 0.23. In both the cases, the best fitted m values are positive, 
but consistent with within la uncertainty, while in the latter case the estimated value for 7 of 1.35 is consistent 
within the errors with the canonical value for radiation. Within the uncertainties, the derived zt are consistent with 
the typical values in literature (e.g. [10, 32]), although, on average, higher. In Fig. 1 we plot the T — z diagram with 
the different best fit curves. 

Assuming a constant value for the ratio '^ /inH — (3, with the constraint < /3 < 1, Lima et al. [16] found as a 
result instead of our eq.(22) the relation 

T{z) = To{l + zy-^ . (32) 

For a value of /3 different from zero the temperature of the expanding universe at high values of z is slightly lower 
than in the standard photon-conserved scenario. Lima et al. [16] discussed the various upper limits which can be 
derived on T{z). The best values are provided by absorption lines from molecules. The results are presented in Fig. 1 
of their paper, from which one can see that the upper limits are well compatible with the standard photon-conserved 
scenario, i.e. (3 = Q. Luzzi et al. [29] found using their data (as given in our Table 1) and also the values mentioned 
in eqs. (23)-(26) as best fit the value 13 — 0.024J^qq24. We performed, as a check, the same fit using the same data 
with in addition the value given in eq.(27) and got (3 — 0.027lo'o27i which clearly is consistent with their value. 
Fitting several data sets to constrain the H{z) function (as given by eq. 19), Ma [15] obtained m = — 0.09lg'°^ and 
^mo — 0.29lg'QY. This model with negative m is plotted as well in Fig. 1 as a comparison, although it is ruled out 
by our fitting of the CMB. 

V. CONCLUSIONS 

We studied a model for the dark energy decay based on the assumption that the dark energy is regarded as a second 
fluid component transferring energy continuously to the material component as suggested by Lima [1] . In which case 
the second law of thermodynamics constrains the whole process, thus allowing to determine the temperature redshift 
relation. As a new point we were able to derive for the latter quantity a relation depending on parameters such as 
Weff and 7. Thanks to recent CMB temperature measurements via the Sunyaev-Zcl'dovich effect and others at higher 
redshift from quasar absorption lines, we could determine best fit values for these parameters. In particular, assuming 
the canonical 4/3 value for 7 we found m = 0.09 ± 0.10 corresponding to Weff = —0.97 ± 0.034, whereas leaving 7 
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FIG. 1: Cosmic microwave background temperature for different m-models. The data are plotted as black symbols. As shown 
in the legend, the gray line and the lighter gray region are the best fitted curve and la scatter for the case with 7 — 4/3 and m 
free. The black dashed line with the cyan region are relative to the case with both 7 and m free. Finally the red curve is the 
best fit curve when using the value for m as found by Ma [15]. In the inset, we show the difference between the data points 
and each best fitted curve as a function of redshift. 

open we got m = 0.20 ± 0.10 or Weff = —0.93 ± 0.08. Both results are within la consistent with the canonical value 
of —1 for the cosmological constant. Using SNe la data combined with CMB and Baryon Acoustic Oscillations, and 
assuming a flat universe, Kowalski et al. [33] found Weff = —0.97 ± 0.06 ± 0.06 (stat, sys) and similarly Kessler et 
al. [34] estimated We/f = —0.96 ± 0.06 ± 0.12. Our best fit value for w^ff compares quite well with these values, and 
we point out that our "derivation" is completely different and to our knowledge new. The values zy we find for the 
transition redshift are in good agreement within the uncertainties with the typical ones quoted in the literature. 

Clearly, it would be nice to have better measurements for the CMB temperature at high redshift, since the difference 
between the different models increases for high z values. We notice also that since the best fit value for m is small, it 
implies that the dark energy variation with time is very slow (at least when not considering too high z values, where 
we expect our approximations to break down) and thus our assumption on w^ — —1 can be considered as justified a 
posteriori. 
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